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ABSTRACT 



We reproduce the blue and red sequences in the observed joint distribution of 
colour and magnitude for galaxies at low and high redshifts using hybrid N-body/semi- 
analytic simulations of galaxy formation. The match of model and data is achieved 
| by mimicking the effects of cold flows versus shock heating coupled to feedback from 

active galactic nuclei (AGNs), as predicted by Dekel & Birnboim (2006). After a 
critical epoch z ~ 3, only haloes below a critical shock-heating mass M s bock~ 1O 12 M 
\Q , enjoy gas supply by cold flows and form stars, while cooling and star formation are 

shut down abruptly above this mass. The shock- heated gas is kept hot because being 
dilute it is vulnerable to feedback from energetic sources such as AGNs in their self- 
regulated mode. The shutdown explains in detail the bright-end truncation of the 
blue sequence at ~ L* , the appearance of luminous red-and-dead galaxies on the red 
sequence starting already at z ~ 2, the colour bimodality, its strong dependence on 
environment density and its correlations with morphology and other galaxy properties. 
Before z ~ 2 — 3, even haloes above the shock-heating mass form stars by cold streams 
penetrating through the hot gas. This explains the bright star- forming galaxies at 
z ~ 3 — 4, the early appearance of massive galaxies on the red sequence, the high 
cosmological star-formation rate at high redshifts and the subsequent low rate at low 
| redshifts. 
& ■ 
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1 INTRODUCTION 

Following the early indications for correlations between 
global galaxy properties alo ng the Hubble sequence (e.g. 
lHubblell92d : lHumasonll936l) . large statistical surveys, such 
as the Sloan Digital Sky Survey (SDSS) and the Two De- 
gree Field (2dF), have revealed the details of a robust bi- 
modality, which divides the galaxy population into a "red 
sequence" and a "blue seque nce". The bimodality is seen in 
(a) the luminosity function (Baldrv et al J 12004 IBell et alJ 
2003), (b) the colour distribution and t he joint distribution 
of colour and magnitude or stellar m ass llBaldrv et alj|2004t 
iBalogh et alJl2004l : iHogg et all2004). (c) the st ellar age and 
the star formation rate iKauffmann et al. 200J3), (d) the gas- 
to-stellar mass ratio dKannappanl 12004) . (e) the bulge-to- 
disc ratio iKauffmann et aTi2003TT (f) the environment den- 



sity telanton^t and (g) the presence of an X-ray 

halo jHelsdon fc PonmaJ2003l : iMathews fc Brighentill2003l : 
lOsmond fc PonmarJl200l) . 

In the colour-magnitude diagram, the red population 
the and blue population form two stripes. In both sequences, 
redder galaxies tend to be brighter, but more so for the blue 
sequence, which is also broader. The blue sequence is trun- 
cated at Al r ~ —22, while the red one extends to brighter 
magnitudes. The division between the two classes of galaxies 
is associated with a critical stellar mass M, crit ~ 3xlO lo M . 
Galaxies below M* C rit are typically blue, star-forming spi- 
rals and reside in the field. Galaxies above M* cr it are domi- 
nated by red spheroids of old stars and live in dense environ- 
ments. The faint r ed galaxies in dense environments define a 
third population felanton et aljl200tJ) . This indicates that 
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the mass of the host halo, rather than the stellar mass of 
the galaxy, determines whether a galaxy is blue or red. 

Already at z ^ 1, the most luminous galax ies form a 
red population separated fro m the blue sequence feell et alJ 
120041 : iMoustakas et al.ll2004l) . whereas at z ~ 2 - 4 there is 
intense star formation in surprisingly big galaxies (Ly man- 
break galaxies and SCUBA sourc es : Sha pley et aD l2004t 
ISmail et alJl2002l IChapman et alJ 120031 l2004h . These ob- 
servations indicate that star formation in massive galaxies 
happened relatively early and then shut off. 

Current models of galaxy formation have difficulties in 
explaining the bimodality features seen at the bright end. 
There is no natural explanation for the low rate of gas supply 
indic ated by the low star form ation rates j Kaiffim ann et all 
2003) and the low HI fractions llKannappanl2 004') of massive 
red galaxies compared to the the high accretion rates and 
gas fractions inferred for blue spiral galaxies. 

In the traditiona l scenario of galaxy formation 
jRees fc Qstrikerl Il977t Isilkl Il977t IWhite fc Reel Il978t 
iBlumenthal et al.lll984l IWhite fc Frenklll99ll) . the infalling 
gas is shock-heated to the virial temperature of the dark halo 
near the virial radius. The inner gas that can cool radiatively 
on a free fall time scale falls in and forms stars, while the 
outer gas is held in quasi-static equilibrium until it eventu- 
ally cools. The "cooling radius" equals the virial radius to- 
day for haloes of mass Af coo i ~ 3x 10 13 (Z/Z Q )/ b M Q , where 
Z is the hot gas metallicity and /b is the cosmic baryon frac- 
tion. In more massive halos the supply of cold gas gradually 
stops. 

This argument has been successful in explaining the 
ballpark of the mass scale separating galaxies from groups, 
but when it is applied to compute galaxy colours, it predicts 
that massive E and SO galaxies should lie on a continuous 
extension of the blue sequence rather than forming a sepa- 
rate sequence that is systematically redder by u — r ~ 0.2 at 
M r H — 22 as seen in the SDSS. Part of this failure arises 

from the gradual dependence of cooling on mass. 

Analyses in spherical symmetry feirnbqim fc Deke 



2003 , BD03) and cosmological simulations iKeres^t^rt 



2005, K05) have shown that in haloes below a critical shock- 
heating mass M s hock ~ 10 12 Mq rapid cooling does not allow 
gas compression to support a stable virial shock. So the gas 
flows cold and unperturbed into the inner halo. The actual 
critical M 8 hock is smaller than the crude cooling mass by a 
factor of a few, and is comparable to the halo mass indi- 
cated by the observed bimodality. The new picture near the 
critical scale is of a multiphase intergalactic medium (IGM) 
with a cold phase embedded in a hot medium. 

These studies indicate that the transition from efficient 
build-up by cold flows to complete dominance of the hot 
mode occurs relatively shar ply o ver a narrow mass range 
about Mshock. iBmnevI J2004h and lDekel fc Birnboiml fcOOrl 
DB06) noticed that this transition can be sharpened further 
once shock-heating triggers another heating meachanism 
that compensates for subsequent radiative losses. Feed- 



back from active galactic nuclei 
candidate for the source of such 


(AGNs) is the leading 


lieating ^il^o^fe Binncv 


19931: iTucker & Davidl Il997t 


Ciotti fc Qstrikerl Il997t 


Silk & Recs 1998: iGranato et al. 


2004: Ruszkowski et alJ 


20041: iBriieeen et alJ 1200,4 iDi Matteo et alJ 1200,4 and 

references therein). The role of AGN feedback is sup- 



in early-type galaxies (7'i lTremaine et all 120021) . (b) the 
evidence from X-ray data that outflows can create deep 
cavit i es in the hot IGM of galaxy clusters jFabian et alJ 
120031 iMcNamara et afll2005l) . (c) the estimate that a small 
fraction of the total energy radiated by an AGN over a 
cosmological time is sufficient, and (d) the possibility of 
a radiatively inefficient AGN mode at low accretion rates 
that allows pumping energy into the IGM without violating 
the strong constr aints on the amount of blue light that can 
be released (e.g. IDi Matteo et al]|2003l . for the particular 
case of M87). 

The critical halo mass for effective feedback does not 
seem to have its natural roots in the physics of the AGNs 
themselves. On the other hand, this scale arises naturally as 
the critical scale for shock heating. Near the shock-heating 
scale, the IGM is made of cold dense clumps in a hot dilute 
medium. This dilute medium is vulnerable to AGN feedback. 
The abrupt shutdown in the accretion of gas by galaxies 
above Mshock is trigerred by shock heating and maintained 
over time by AGN feedback. 

DB06 have studied the potential implications of the in- 
terplay between the cold/hot flow transition and different 
feedback sources and have predicted how this can reproduce 
the observed bimodality features. In this paper, we test this 
hypothesis in a more quantitative way, by incorporating the 
shutdown above a critical mass into simulations of galaxy 
formation that combine N-body gravity with semi-analytic 
modelling (SAM) of the baryonic processes. In i|5|we present 
the semi-analytic/N-body code. In |H]we show that without 
the abrupt shutdown of gas accretion above a critical mass 
the SAM overpredicts the number of bright blue galaxies 
and fails to reproduce the red colours of massive early-type 
galaxies. In f|Hwe summarise the theory of cold flows versus 
shock heating, and describe how we implement these two 
features in the SAM. In we explain our picture of the 
interaction between the massive black hole and the IGM. In 
we find that the addition of such an abrupt shutdown 
brings the model into excellent agreement with the data. In 
33 we demonstrate that the success of this model predom- 
inantly depends on the requirement that the shock-heated 
gas never cools and is robust to the detailed implementation 
of other physical ingredients in the model. In fJUwe interpret 
the results of our simulations and propose a unified scenario 
for the origin of the blue sequence and the two parts of the 
red sequence. In J3|we summarise our main conclusions. 



ported by (a) the ubiquity of supermassive black holes 



2 THE HYBRID SEMI-ANALYTIC/N-BODY 
CODE 

GallCS (Galaxie s In Cosmological Simulations; 
lHatton et alJ [2003 ) is a method that combines high 
resolution simulations of gravitational clustering of the 
dark matter with a semi-analytic approach to the physics of 
the baryons (gas accretion, galaxy mergers, star formation 
and feedback) to simulate galaxy formation in a ACDM 
Universe. The version of GallCS used for this article differs 
from the original version in the star formation law (Eq. 1), 
the outflow rate due to supernova feedback (Eq. 2) and the 
cut-off that prevents cooling on massive galaxies (Eq. 3). 
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2.1 The dark matter simulation 

The cosmological N-body simulation has been carried out 
with the parallel tree code developed by Ninin (1999). The 
cosmological model is flat ACDM with a cosmological con- 
stant of £Ia = 0.667, a Hubble constant of Ho — 66.7 km s -1 , 
and a power spectrum normalised to as = 0.88. The simu- 
lated volume is a cube of (150 Mpc) 3 with 256 3 particles of 
8.3 x 1O 9 M0 each and a smoothing length of 29.3 kpc. The 
simulation produced 100 snapshots spaced logarithmically 
in the expansion factor (1 + z) -1 from z = 35.59 to z = 0. 

In each snapshot, w e applied a friend-of-friend algo- 
rithm jDavis et alj Il985l) to identify virialised haloes of 
more than 20 particles. The minimum halo mass is thus 
1.65 x 1O 11 M0. The information about individual haloes 
extracted from the N-body simulation and passed to the 
semi- analytic model is in three parameters: the virial mass, 
the virial density and the spin parameter. Merger trees are 
computed by linking haloes identified in each snapshot with 
their progenitors in the one before, including all predeces- 
sors from which a halo has inherited one or more particles. 
We do not use the substructure provided by the N-body 
simulation. Once a halo becomes a subhalo of another halo, 
we switch from following its gravitational dynamics with the 
N-body integrator to an approximate treatment based on a 
semi-analytic prescription. 

2.2 Galaxy formation through the accretion of gas 

Newly identified haloes are attributed a gas mass by as- 
suming a universal baryon fraction of fib/^o = 0.135. All 
baryons start in a hot phase shock-heated to the virial tem- 
perature. Their density profile is that of a singular isother- 
mal sphere truncated at the virial radius. 

The cooling time is calculated from the density dis- 
tribution of the hot gas using the cooling function of 
ISutherland fc Dopital il993f) . The gas for which both the 
cooling time and the free fall time are shorter than the 
timestep At is allowed to cool during that timestep. Cooling 
is accompanied by inflow to maintain the shape of the hot 
gas density profile and is inhibited in haloes with positive 
total energy or with angular momentum parameter A > 0.5. 

A systematic comparison of GallCS and SPH results 
(Cattaneo et al. 2006) shows that the galaxy mass function 
obtained with the GallCS cooling scheme is very similar to 
the one found with the SPH method, which demonstrates 
that the cooling algorithm is not an important source of 
error. 

2.3 Morphologies 

GallCS models a galaxy with three components: a disc, a 
bulge and a starburst. Each galaxy component is made of 
stars and cold gas, while the hot is considered as a compo- 
nent of the halo. Only the disc accretes gas from the sur- 
rounding halo. The bulge grows by mergers and disc insta- 
bilities, while the starburst contains gas in transition from 
the disc to the bulge. The disc is assumed to be exponential, 
with the radius determined by conservati on of angular mo- 
mentum. The bulge is assumed to have a lHernauisd (ll990T) 
profile and its radius is determined based on energy conser- 
vation. 



Disc instabilities transfer matter from the disc to the 
bulge until the bulge mas s is large enough to stabilize the 
disc llvan den Boscbiri998t) . Dynamical friction drives galax- 
ies to the centres of dark matter haloes and is the most 
important cause of mergers (while we also include satellite- 
satellite encounters). The fraction of the disc that trans- 
fers to the bulge in a merger grows with the mass ratio of 
the merging galaxies from zero, for a very minor, merger 
to unity, for an equal mass merger. This is of course a sim- 
plified picture of the dynamics of morphological transforma- 
tions, but it provides results consistent with essential proper- 
ties such as the Faber-Jackson relation and the fundamental 
plane of spheroids. 

2.4 Star formation and feedback 

The star formation law is the same for all components: 

M. = ^211(1 + 0)-. (1) 

The cold gas mass Af co id refers to the component in question, 
and idyn is the dynamical time (corresponding to half a ro- 
tation for a disc and the free fall time for a spheroid) . In the 
standard GallCS model, star formation is activated when 
the gas surface density is E gas > 20m p cm~ 2 (m p is the 
proton mass) . The star for mation efficiency param eter has a 
fiducial value of /3* = 50 iGuiderdoni et alJll998l) . which is 
assumed to be the same at all redshifts, a, — 0. Gas in tran- 
sit from the disc to the spheroid passes through a starburst 
phase in which the star formation rate grows by a factor of 
10 before the gas is converted into bulge stars. This high 
star formation rate is obtained by assuming that the star- 
burst radius is ten times smaller tha n the final bulge radiu s 
(see the hydrodynamic simulation in lCattaneo et alll2005h . 
The properties of the galaxy population at low redshifts are 
insensitive to the starbust star formation time scale as long 
as it is much shorter than the Hubble time. 

The mass of newly for med stars is distributed accord- 
ing to the lKennicutti al983f) initial mass function. Stars are 
evolved between snapshots using substeps of at most 1 Myr. 
During each sub-step, stars release mass and energy into 
the interstellar medium. Most of the mass comes from the 
red giant and the asymptotic giant branches of stellar evo- 
lution, while most of the energy comes from shocks due to 
supernova explosions. The enriched material released in the 
late stages of stellar evolution is mixed with the cold phase, 
while the energy released from supernovae is used to reheat 
the cold gas and return it to the hot phase in the halo. Re- 
heated gas is ejected from the halo if the potential is shal- 
low enough. The rate of mass loss through supernova-driven 
winds M w is determined by the equation 

l- 2 

-M wj u csc = csn?7sn-EsnM», (2) 

where _Esn = 10 51 erg is the energy of a supernova, t]sn = 
0.0093 is the number of sup ernovae every 1 Mp , of stars and 
Vcsc is the escape velocity ijDekel fc Silklll98m . In GallCS 
v csc ~ 1.84i> c for discs and w esc = 2a for bulges/starbursts. 
The supernova efficiency esn is a free parameter, which de- 
termines the fraction of the supernova energy used to reheat 
the cold gas and drive an outflow at the escape speed. We 
determine its value by requiring that the model matches the 
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luminosity and the gas fraction of a galaxy that lives in a 
Milky Way-type halo. This constraint is consistent with the 
observed normalization of the Tully-Fisher relation (Fig.0. 
Our best fit value esn — 0-2 agrees with those adopted i n 
other SAMs JSomerville fc Primacldll999HCole et alJl200Cft . 

2.5 Stardust 

The STARDUST model foevriendt et al.l ll99Sh gives the 
spectrum of a stellar population as a function of age and 
metallicity and includes a phenomenological treatment of 
the reprocessing of light by dust. GallCS uses STARDUST 
to compute the spectrum of each component and output 
galaxy magnitudes. Dust absorption is computed with a phe- 
nomenological extinction law that depends on the column 
density of neutral hydrogen, the line of sight and the metal- 
licity of the obscuring material. The reemitted spectrum is 
the sum of four templates (big and small carbon grains, 
silicates, and polycyclic aromatic hydrocarbons). Their rel- 
ative weights are chosen to reproduce the relation between 
bolometric luminosity and infrared colours observed locally 
in IRAS galaxies. In GallCS we consider only the self- 
absorption of each component. The mean column density is 
determined assuming spherical symmetry for the starburst 
and the bulge, while in the disc component there is a depen- 
dence on a randomly selected inclination angle. 



3 FAILURE OF THE STANDARD MODEL 

Semi-analytic models based on the standard scenario (SJ5J 
have difficulties in explaining the early formation of bright 
elliptical galaxies followed by passive evolution. At z = 0, 
the luminosity function shown in Fig. [5] demonstrates that 
the standard model predicts too many bright galaxies com- 
pared to SDSS. 

Figure[3]and Fig. 0] illustrate in two different ways that 
the model bright ellipticals are also not red enough. In Fig. [3] 
we see that the brightest galaxies (which are also the central 
galaxies of massive haloes) fall below the colour range of red 
galaxies in the SDSS. Figure 2] displays the galaxy colour 
distribution in different magnitude bin. A comparison of the 
predictions of the standard model with the SDSS data in the 
—23 < M r < —22.5 bin shows that (a) there are many more 
galaxies in the model than in the data, (b) the predicted 
peak of the colour distribution is bluer by u — r ~ 0.2, and 
(c) the model distribution is skewed toward blue colours. 

At z ~ 1, Fig. |S] shows that the model fails to predict 
the existence of bright red galaxies. Instead, the most lumi- 
nous galaxies are predicted to lie on an extension of the blue 
seque nce, in conflict with the COMBO- 17 data (Be ll et all 
2004), which reveal the bimodality already at z ~ 1. 

Fig. [U shows the failure of the standard model in 
mat ching the luminosi ty function of Lyman-break galax- 
ies JSteidel et alJll999h at z ~ 3. With its star formation 
rate, the model cannot explain the presence of massive star- 
forming galaxies at z ~ 3. In the whole computational vol- 
ume, there is only a handful of galaxies with Mv < —23 
(Fig- 13 ■ The star-formation history shown in Fig. [^demon- 
strates again that the model star formation is not sufficient 
at high redshifts, and that it is too high at low redshifts. 
These discrepancies are generic to the standard scenario and 



are not specific to the way this scenario is implemented in 
GallCS. iBlaizot et all ll2004h were able to use GallCS with- 
out the modifications propo sed here and they di d find a rea- 
sonable agreement with the lSteidel et al] lll999F) luminosity 
function at z ~ 3, but they used a non-standard feedback 
model, which spoils the fit to joint distribution of colour and 
magnitude at z ~ 0. 



4 COLD FLOWS VERSUS SHOCK HEATING 

A stable extended shock can exist in a spherical system when 
the pressure that develops in the post-shock gas is sufficient 
to balance the gravitational attraction toward the halo cen- 
tre. For adiabatic contraction, this requirement translates 
into a condition on the adiabatic index 7 ~ 5/3, defined as 
the ratio between the specific heats at constant pressure and 
constant volume. 

BD03 have extended this result to the case where energy 
is lost at a rate q by replacing the adiabatic index 7 with the 
effective polytropic index 7 c fr = dlnp/dlnp = 7— (p/ p)(q/e), 
where e is the internal energy. This can be rewritten as jgg = 
7 — icomp/tcooi, where t coo \ = e/q is the cooling time of the 
post-shock gas and t CO mp = p/p = — (Vu) -1 = — r s /(3ui) is 
its compression time, with r a the shock radius and u\ the 
post-shock radial velocity. The first equality is the continu- 
ity equation and the second assumes that the post-shock 
radial flow pattern is homologous. Using perturbation anal- 
ysis, BD03 showed that the criterion for a stable shock is 
7 c ff > 10/7 (for a monoatomic gas in spherical symmetry), 
which translates to t~^ ol < t" 1 ^ . This defines a critical halo 
mass, above which the gas is shock-heated near the virial 
radius, and below which the gas flows cold and unperturbed 
into the inner halo, where it may eventually shock. This 
result has been confirmed by BD03 using spherical hydro- 
dynamical simulations. A similar conclusion was obtained 
earlier i n one- dimensional simulations of pancake formation 
l|Binnevlll977h . 

DB06 have applied this theory in a cosmological con- 
text. They have assumed that the pre-shock gas traces the 
dark matter distribution in its free fall, and the strong 
shock condition determines the post-shock infall speed, den- 
sity and temperature. They infer a shock-heating scale of 
A/shock ~ 6 x 10 11 Mq for a shock in the inner halo (r ~ 
0.1r\,i r ), and Af s hock ~ 3 x 1O 12 M0 for a shock at the virial 
radius. 

Cosmolog ical simulations w ith smoothed particle hy- 
drodynamics IIFardal et alJl200lt K05) and with an adaptive 
Eulerian mesh (DB06) have shown that the phenomenon is 
not an artifact of the spherical or planar symmetry. Fig. 6 
of K05, based on a simulation with zero metallicity, shows 
that the infall is almost all cold at M^aio ^ 10 M@, drops 
to 50% at M s hock ~ 3 x 10 U M Q , and for z < 2 is almost 
all hot at M ha i > 10 12 M o . This is an underestimate of 
the critical mass by a factor of order two because of the as- 
sumed zero metallicity. At higher redshifts the cold fraction 
does not vanish even above the critical mass. At z = 3 it 
is ~ 40% even in the most massive halo contained in the 
computational box, 1O 13 M0. 

We model the complex multiphase behavior in an ide- 
alised way: a sharp transition from 100% cold flow to a 
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Standard Model New Model 




11 12 13 14 15 11 12 13 14 15 

M ho ,„ (M Q ) M holo (Mo) 

Figure 1. Gala xy luminosity versus halo mass for the "standard" model (left) and the "new" model (right). The blue square refers to 
the Milky Way feinnev fc Merr ificld 199 & IPehnen &; Binne\lll99g|) . The vertical dashed line marks the shock-heating scale. The main 
difference between the two models is in the break in the relation at M ha j > M shock due to the shutdown above M cr ; t (Eq. 3). Some 
differences also show up below Afghocls because of the shutdown of cooling in galaxies where the bulge is the dominant component. 



Standard Model New Model 




-24 -23 -22 -21 -20 -24 -23 -22 -21 -20 



Figure 2. Luminosity function in the r— band. The model predictions (solid line, red) for the "standard" model (left) and for the "new" 
mode l with shutdown in massive haloes (right) compared to the observed luminosity function (symbols) observed by SDSS ( Baldrv et al. 
|2004|) . The dotted-dashed line (right) illustrates the effect of lowering the critical redshift from z c = 3.2 to z c = 3. 




Figure 3. Colour— magnitude diagram at 2 = using u — r versus M r (SDSS magnitudes). Mode l predictions ( dots), "standard" (left) 
and "new" (right), compared to the observed loci of the red/blue sequence in the SDSS data fealdrv et all2004l) marked by the pairs of 
red/blue curves. Galaxies in haloes below Mg^ock are marked blue; they are typically the only galaxy in their halo. Galaxies in haloes 
above M ghn r |<- are marked red if they are the central galaxy of their halo and green if they are satellites. The contours, directly comparable 
to Fig. 2 of lBaldrv et alj (2004)., refer to the density of points in units of mag 2 Mpc 3 ) and highlight the bimodality. 
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Figure 4. u — r colour distribution in r-band magnitude bins. The model predicti ons (red histogram), "standard" (left) and "new" (right), 
compared to the SDSS data (black smoothed histogram) from lBaldrv et al.U2004l) , The "new" model parameters are M B[loc ^ = 2xlO 12 M0 
and z c = 3.2 (eq. |3J. The star formation efficiency increases at high redshift with a power of a„ = 0.6 (eq. and the energetic efficiency 
of supernova feedback is egN = 0.2 (eq. In- 



complete shutdown of the accretion at a critical halo mass 
M CT it {z). We parameterise the critical mass by 

M crit = Af shock x min{l, 1 o 1 - 3 («-«<>>} . ( 3 ) 

At redshifts below z c the critical mass is the shock-heating 
mass, which is expected to be of the order of 1O 12 M0. We 
treat its exact value as a free parameter determined by the 
best fit to different constraints. The predicted presence of 
cold flows above the shock-heating mass at z > z c is mim- 
icked by the exponential increase in Merit- The critical red- 
shift z c , predicted in the ballpark of ~ 2 — 3, is also adjusted 
for best fit. 



5 NEW ELEMENTS IN THE GALAXY 
FORMATION SCENARIO 

5.1 shutdown by virial shock heating and AGN 
feedback 

The critical stellar mass of 3 x 1O 1O M0 associated with the 
observed bimodality corresponds to a critical halo mass of 



10 Mq. Star forming galaxies are mainly confined to 
haloes below this mass and most galaxie s in haloes above 
this threshold are evolving passively (Ka uffmann et all 
2004). The main reason for the failure of the standard model 
at low redshift is the cooling of hot gas onto the central 
galaxies of haloes that are more m assive than the critica l 
scale. The traditional cooling scale l|Rees fc Ostrikerlll97'7l) 
is somewhat too big and is not associated with a transition 
that is sharp enough to explain the sharp colour bimodality 
and the steep drop of the luminosity function beyond L* . 

The gas that cools and accumulates in the central 
galaxies of massive haloes will inevitably form stars un- 
less there is a mechanism that reheats it or removes 
it. AGNs p rovide the most plau s ible source of energy 
for the job jBinnev fc TaboJll995t lTucker fc Davidl 119971: 
Ciotti fc Ostrikerlll997USilk fc Reejll 994 iNnlsen fc Fahia.nl 
2000; iSpringel et alJl2005T) . The supporting evidence comes 
from several lines of argument such as (a) the abun- 
dance of supermassive black holes in massive spheroids to- 
gether with the correlation between black hole mass and 
bulge mass jMarconi fc Huntll2003l lHaring fc Rb3l2004l) or 
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the bulge's velocity d ispersion jMerritt fc Ferraresd 1200 ll 
iTremaine et alJ l2002h . (b) the concordance between the 
quasar epoch and the stellar ages of early-type galaxies 
jGranato et alJ l200ll : ICattaneo fc Bernard! 120031) . (c) the 
connectio n of black hole growth and elliptical gala xies with 
mergers jToomre fc Toomrdll972l IStocktonlll999D . and (d) 
the observation of powerful outflows from active galaxies 
Ce.g.. lMcNamara et aT1l2005l) . 

Supernova feedback is not powerful enough to affect the 
gas in massive haloes significantly unless the rate of for- 
mation of massive stars is especially high, but then these 
galaxies would no longer be red. A similar argument can ap- 
ply to feedback from quasars since since early-type galaxies 
are 10 4 times more c ommon than quasars at low redshifts 
iWisotzki et al.ll200lfl . This is a serious problem for mod- 
els in which AGNs are supposed to interact with the halo 
gas through radiative processes such as radiation pressure 
or Compton scattering. However, there is evidence show- 
ing that AGN outflows are possible not only during the 
main phase of black hole growth, which was very brief and 
happened at high redshift, but also when the black hole 
accretion rate is so low that the AGN is not optically lu- 
minous (see e.g. the adiaba tic inflow-outflow solution by 
iBlandford fc Begermar]|l999l and the study of the jet in M87 
bv lDi Matteo et al.ll2003F ~ In the inflow-outflow model, the 
density of the accretion flow is too low for the gas to radi- 
ate efficiently. The accretion power is released mechanically, 
through particle jets, rather than radiatively, through the 
emission of optical/UV light. In this scenario, one possible 
way for the AGN to heat the gas is through shocks pro- 
duced by the propagation of jets (e.g. iRevnolds et aljEoOll 
iFabian et ai1l2003l: lOmma et al.ll2004h . although the prob- 
lem of how the heat generated by this pr ocess is distributed 
in the hot gas remains open (see also iBegelman fc Nathl 
120051 : IFabian et alll2005l and references therein). 

One difficulty in connecting AGN feedback to the 
galaxy bimodality is that the critical scale does not seem 
to affect black hole growth. For example, the correlation be- 
tween black holes and bulges extends down to small bulges 
as in the Milky Way (iMarconi fc Hundl2003l; lHaring fc Rb] 
I2004I) . while cooling seems to shut down only in massive 
galaxies. 

Our proposal here (and in DB06) is that AGN feed- 
back is switched on by the change in the large scale black- 
hole environment occuring once the gas is shock heated 
above the shock-heating scale. While the cold phase is frag- 
mented in dense clouds, the hot phase is distributed much 
more uniformly in a low density medium. Outflows from 
the AGN will preferentially propagate through the less re- 
sistive hot gas that fills the space between the clouds and 
leave the cold clouds behind instead of blowing them away 
(e.g.. lMcKee fc Ostrikeill977l. in the case of supernova blast 
waves.). The dilute hot gas is thus more vulnerable to the 
propagation of shock fronts. This serves as the basis for the 
new scenario simulated in the current paper. 

We thus adopt the assumption that all massive galax- 
ies contain a supermassive black hole, which can accrete gas 
and deposit energy in the surrounding gas, but the capacity 
of the gas to react to this injection depends on the presence 
of a shock-heated component. As soon as this condition is 
fulfilled, the gas begins to expand and the black hole accre- 
tion rate drops until it stabilises at a value that over time 



compensates the radiative losses of the hot gas. The onset of 
this self-regulated accretion cycle prevents the shock-heated 
gas from ever cooling. We model this physical scenario by 
interpreting the critical mass that separates the regimes of 
cold and hot accretion (eq. |2J as a sharp cooling cut-off. 

When most of the mass is in the hot phase, the cold 
clouds also become vulnerable to feedback. We therefore 
assume that when Mhaio > M cr i t the cold gas in the cen- 
tral galaxy is reheated to the virial temperature of the halo. 
Black hole and supernova winds from gas-rich galaxy merg- 
ers will also sharpen the transition to a hot IGM because the 
cold streams in free fall will go through a shock when they 
hit the winds blown from the galaxy. We account for this 
scenario by preventing cold flows on galaxies in which the 
bulge is the main component as a massive bulge correlates 
with a massive black hole and, therefore, with a history of 
high AGN activity. This bulge-related condition improves 
the details of the quantitative agreement with the SDSS 
colour-magnitude distribution, but we will see below (SjTJl 
that it plays only a minor role in curing the discrepancies 
and achieving a good match with the observed features. 

5.2 Cold streams and star formation at high 
redshifts 

The standard scenario under-predicts the number of bright 
galaxies in the rest-frame 0.17 pm luminosity function at 
z ~ 3 (Fig. This is due to a slow conversion of gas into 
stars, as it is not cured even when all the gas is allowed to 
cool in free fall from rest and supernova feedback is practi- 
cally shut off. The shutdown of cooling in the most massive 
haloes will only accentuate the problem. The problem can 
be cured if the star formation recipes derived from local ob- 
servations are revised to take into account the different way 
by which galaxies accrete their gas at high z. As discussed 
in 21 t ne cold gas is supplied by streams at the virial ve- 
locity along the dark matter filaments. The collision of the 
cold streams among themselves and with the galactic disc 
produces bursts of star formation analogous to the bursts 
resulting from the collision of two disc galaxies or cold gas 
clouds. Under the conditions that allow a cold flow, these 
collisions are expected to produce isothermal shocks, and 
the rapid cooling behind the shocks generates dense cold 
slabs in which the Jeans mass becomes small and stars can 
form efficiently. 

The detailed physics of star formation under these con- 
ditions is yet to be worked out, but this argument strongly 
suggests that the conversion of gas into stars must have been 
more efficient at high z than it is now. We mimic this effect 
in the "new" model by introducing a ~ (1 + z) a * increase 
in the star formation efficiency (eq. 0, where a* is a free 
parameter determined by fitting the luminosity function of 
Lyman-break galaxies at z ~ 3. 



6 RESULTS 

We have revised the GallCS SAM to accommodate the new 
features discussed above, in particular (a) the high star for- 
mation rate due to cold streams at high redshift, (b) the 
abrupt shutdown of galaxy accretion and star formation due 
to shock heating in haloes above a critical mass, and (c) 
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Figure 5. Colour-magnitude diagram at z = 1. U — V versus M^. The "standard" and "new" models are the same as in Fig.0 The 
red solid line and the black dashed line refer to the centre and the lower bound of the red sequence for the 1.0 < z < 1.1 bin in the 
COMBO-17 data feell et al.ll2004h . The dot-dashed line marks the estimated completeness limit of that survey. 
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Figure 6. Luminosity function at z ~ 3.1. The magnitude is rest-frame 0.17 /im. The model predictions (solid, red line), "standard" (left ) 
and "new" (right), compared to the data (symbols) inferred from ij-band observations of Lyman-break galaxies JSteidel et al.lfl999tl . 
Superimposed are data from the VIMOS/VLT Deep Survey (VVDS). The dashed line in the left panel refers to an extreme variant of 
the "standard" model where all the gas cools instantly and feedback is practically shut off. Even with this extreme cooling the model 
fails to form enough bright star-forming galaxies at z ~ 3, demonstrating that a higher star formation efficiency is necessary at high z. 
The dot-dashed line in the right panel refers to the "new" model with z c = 3.0 rather than the fiducial z c = 3.2, showing that the high-z 
predictions of the "new" model are sensitive to the critical redshift after which M cr i t = iW s h ck- The thin dotted line in the right panel 
corresponds to 2 C ^ oo, namely M CI n = M s i! OC k for all z. 



the maintenance of this shutdown by AGN feedback that 
couples to the hot gas. This involves three new free SAM 
parameters, which are fine-tuned to reproduce the obser- 
vational constraints. They are: (a) the shock-heating mass 
scale, Mshock, (b) the critical redshift z c prior to which M cr i t 
is growing with z (eq. |HJ , and (c) the power index a* in 
the redshift dependence of the star formation rate (eq. 0. 
The shock-heating scale is the single most important pa- 
rameter determining the fit at low and intermediate red- 
shifts. As seen in the M r -Mh a io diagram of Fig. the crit- 
ical halos mass translates into a maximum luminosity for a 
star-forming galaxy, where the main inclined stripe mark- 
ing the brightest galaxies intersects the vertical red line. 
This parameter is tuned to match the number density of 
blue (u — t £ 2) galaxies with M r ~ —22.25 as observed 
in the SDSS. The other two parameters are determined by 
fitting the luminosity function of Lyman-break galaxies at 



dSteidel et alJll999^ . The obtained best fit values are 
M shock = 2 x 10 12 M Q , z c = 3.2, a, = 0.6. 

With this choice of the parameter values, we recover an 
excellent match to the z = luminosity function (Fig.^J and 
to the bimodal colour-magnitue distribution in the SDSS 
(Fig. E1& Fig. 2J. Figure El shows that the blue sequence is 
properly truncated at M r ~ —22.25. The central galaxies 
of massive haloes are also positioned correctly in the SDSS 
red sequence colour range. The colour histograms (Fig. 
show that the simulated joint colour-magnitude distribu- 
tion agrees with the observed distribution both in shape 
and in normalisation. We now reproduce the red population 
observed at z ~ 1 in the COMBO-17 survey and the fact 
that the brightest galaxies are indeed on the red sequence 
(Fig. 0. Finally, this model provides an excllent fit to the 
comoving number density of bright star-forming galaxies at 
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Figure 7. Colour-magnitude diagram at z = 3. The "standard" and "new" models are the same as in Fig. The lines, shown for 
reference only, are the same as in Fig. ^describing the COMBO-17 red sequence at z ~ 1. 
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Figure 8. Global star-formation history. The model pred ictions (solid, red line), "standard" (left) and "new" (right), compared to the 
data (symbols) from the GOODS survey iGiavalisco et all20o3) . While the "standard" model overpredicts the star-formation density at 
z < 2 the "new" model provides an excllent match to the data at all redshifts. 



z ~ 3 (Fig. [SJ and to the observed cosmological star forma- 
tion history (Fig. |HJ . 



7 SENSITIVITY TO MODEL INGREDIENTS 

Here we explore the contribution of different ingredients 
to the success of the model in fitting the observed colour- 
magnitude distribution. 

7.1 The critical mass at low z 

The shutdown above the critical halo mass is the main 
new feature of our model, and thus M srl ock is the key pa- 
rameter. Figure compares the distribution of galaxy lu- 
minosity versus halo mass in the "standard" model and 
our "new" model. Preventing cooling and star formation in 
haloes above M s h oc k keeps the magnitudes of these galaxies 
fainter than M r ~ — 23 The abrupt change in the mass-to- 
light ratio at the critical scale separates the red sequence 
from the blue sequence and sets an upper limit to the lumi- 
nosity of galaxies in the blue sequence. The best fit to the 



data, with M shock = 2 x 10 12 M Q , is shown in Fig. 21 Fig- 
urc[|5]shows how by lowering (or raising) Mshock to 1O 12 M0 
(or 3 x 10 12 Mq) the model reproduces too few (or too many) 
galaxies at M r — —22.25. 



7.2 The critical mass at high z 

At high z, our fiducial "new" model with z c ~ 3 allows star 
formation even in haloes above M s hock (eq. |3J , mimicking 
cold streams in hot haloes. This change is responsible for 
the appearance in Fig. [7| of a population of bright galax- 
ies with a high star forma tion rate at z ~ 3 (compare to 
Sawi cki fc Thompson! l2005l while the red sequence of pas- 
sively evolving galaxies has not formed yet. 

The model predictions for Lyman-break galaxies 
(LBGs) are sensitive to the value of z c . When the increase 
of Merit at high redshifts is not taken into account, the pre- 
dicted luminosity function at 2 ~ 3 falls short at the bright 
end compared to that of bright LBGs (z c — > 00 in Fig. |SJ. 
On the other hand, when z c is taken to be smaller than 
the fiducial value of 3.2, the predicted bright end of the lu- 
minosity function overshoots that of LBGs as observed by 
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ISteidel et alJ (Il999l) . If, however, the actual comoving den- tion at z = becomes even better (Fig. |3J, but in general 
sity of bright LBGs is somewhat higher, as indicated by the the effect of z c on the results at low redshifts is negligible 
VIMOS/VLT Deep Survey (VVDS) team (private commu- (Fig. He and d). 
nication), then z c = 3.0 becomes our best fit value (Fig.|SJ. 
With this value, the match to the r-band luminosity func- 
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Figure 9. The effect of varying different model ingredients on the joint distribution of colour and magnitude at z = 0. Shown is the u — r 
colour distributi on in r-band magnit ude bins for the model predictions (red histogram) compared to the SDSS data (black smoothed 
histogram) from Ifjaldrv et alj J2004T) . The models are the fiducial "new" model of Fig. [3] with one ingredient varied as follows: (a) 
M shock = 10 12M O (compared to 2 X 10 12M © in the fiducial model), (b) A/ ahock = 3 X 10 12 M Q . (c) z c -* do, i.e., same Af crit at all z. (d) 
z c = 3, somewhat lower than the fiducial value, (e) No dominant-bulge requirement for shutdown of cooling and star formation, (f) No 
ejection of the remaining cold gas when cooling is shut off. (g) No increase in star formation efficiency at high rcdshift (a* = 0). (h) An 
increased efficiency of supernova feedback (esn = 0.25). 
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7.3 A bulge criterion for shutdown 

Our fiducial model allows no cooling or star formation in 
galaxies where the bulge is the dominant component. In 
practice, given the resolution limit of the N-body simula- 
tion, the critical shutdown scale is reduced, in this case, to 
the resolution scale of ~ 2 x 10 11 Mq. This scale is signif- 
icantly lower than the best-fit value of M s hock, but is only 
slightly lower than the critical halo mass as predicted by the- 
ory ( BD03; K05; DB06) and as observed jKauffmann et all 
2003). As explained in ^ our M s h oc k is closer to an upper 
limit than to a mean value for the physical shock-heating 
scale. 

When the effect of the bulge fraction is eliminated from 
the criterion for shutdown (Fig. [ij), the model shows a 
small excess of bright galaxies (M r ~ —22.25), both blue 
and red. This is because the more efficient cooling allows 
more massive star-forming galaxies, which in turn speeds 
up the merger rate due to more efficient dynamical friction. 
We learn that the bulge criterion has some effect, but it ba- 
sically serves as a secondary criterion that helps fine-tuning 
the match to the data. The fit is fairly adequate without it. 

7.4 Cold gas ejection in shock-heated haloes 

In the fiducial model, once the halo mass grows above the 
critical mass, all the gas remaining in the central galaxy 
is heated to the virial temperature and all modes of star 
formation are shut off. The idea is that when most of the 
gas is hot, AGN outflows heat, destroy or blow away the 
cold clouds as well. When the cold gas is not ejected from 
massive galaxies but rather allowed to form stars later, the 
model predicts a small excess of blue galaxies in the brightest 
bin (Fig. |Sf ) ■ This is similar to the effect of eliminating the 
bulge criterion, but even less noticeable. 

7.5 Efficient star formation at high redshift 

We have shown that increasing the star formation efficiency 
at high z is necessary to reproduce the luminosity function 
of LBGs. When the star formation is kept constant with red- 
shift, instead, star formation is postponed to later epochs. 
The model predicts a shortage of galaxies in the brightest 
bin and the colours become slightly too blue, but these are 
weak effects (Fig. Ek) . 

7.6 Supernova feedback 

Our fiducial model for supernova feedback is similar to 
the model used in most semi-analytic models based on 
iDekel fc Silkl dl986Tl .lt is more efficient in less massive galax- 
ies with shallower potential wells. The naive expectation is 
that a stronger supernova feedback would make the galax- 
ies redder by removing gas and stopping star formation. 
Instead, we find that supernova feedback makes galaxies 
slightly bluer by ejecting metals from the galaxies into the 
IGM (Fig. Eh) - The metal-enriched gas is later accreted onto 
the hierarchically assembled more massive galaxies, and it 
is kept there as supernova feedback becomes less efficient. 
Thus, supernova feedback ends up transferring metals from 
low mass galaxies into massive galaxies. 



We conclude that the main element responsible for the 
improved match of the model to the colour-magnitude data 
at low redshifts is the abrupt shutdown above a critical mass 
of ~ 10 12 A/q. All the other changes in the model recipes, 
including the bulge criterion, are of secondary importance. 
They serve for fine-tuning the model in order to achieve a 
nearly perfect match. 



8 CORRELATIONS WITH OTHER 
PROPERTIES 

The observed bimodality of the galaxy population, which 
is very apparent in the colour-magnitude diagram, also in- 
volves all other global properties of galaxies. We already 
saw in Fig. [3] that there is a strong correlation between be- 
ing blue or red and being below or above the critical halo 
mass. In Fig. 1101 we explore a variety of such correlations as 
predicted by the "new" model at z — 0. In the same colour- 
magnitude diagram shown in each panel, colour refers to 
a different galaxy property. The properties considered here 
include environment density via halo mass, bulge fraction, 
stellar age, stellar metallicity, total stellar mass and star- 
formation rate. 

A useful way to interpret the various correlations dis- 
played in Fig. 1101 is to consider the evolutionary track of a 
typical galaxy in the colour-magnitude diagram. Panels a 
and e show that the blue sequence is a sequence of growing 
halo mass and stellar mass from bottom- right to top-left. As 
a halo is growing in mass, starting from the smallest mass 
resolved by the N-body simulation, gas is accumulating in 
the disc of the central galaxy and stars are forming. The 
galaxy becomes brighter because its stellar mass increases, 
but it also gets redder because a growing fraction of the 
stellar mass is in old stars. In this evolutionary phase, the 
galaxy moves along a stripe which roughly coincides with 
the blue sequence of the SDSS marked by the blue lines in 
Fig.Eni 

A galaxy reaches the top of the blue sequence when 
its host halo mass becomes comparable to M cr it- Soon af- 
ter, it ceases to accrete gas, star formation shuts off, and 
the galaxy evolves passively to become red and dead. Most 
galaxies never reach the top of the blue sequence because 
they merge into a halo much more massive than their own 
before their original halo exceeds the critical mass. These 
galaxies are typically of low mass compared to the halo into 
which they have merged, and therefore their dynamical fric- 
tion time for sinking into the new halo centre is long. They 
become the "satellite" galaxies populating groups and clus- 
ters of galaxies. Since satellite galaxies are assumed not to 
accrete gas, they fairly quickly exhaust their remaining gas 
and become passive. Their brightness fades and their colour 
becomes redder, as the most massive, bluest stars are the 
first to die. These evolutionary tracks can be deduced, for ex- 
ample, from the stellar mass panel of Fie. 1101 where galaxies 
evolve passively along the equal colour stripes, which stretch 
roughly perpendicular to the blue sequence, from bottom- 
left to top-right. By mentally superimposing these diagonal 
equal-stellar-mass stripes on the panel of Fig. llOl referring to 
halo mass, one can verify that the evolution up along such 
stripes involves an increase in halo mass, consistent with 
merging into a more massive halo. A similar inspection of 
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Figure 10. Correlation of colour— magnitude with other properties in the "new" model at z = 0. As in Fig.|2] the points mark individual 
galaxies in the simulation and the contours refer to their number density in the colour-magnitude plane (mag - 2 Mpc — 3 ), while the blue 
and red lines mark the blue and red sequences in the SDSS. All the above are the same in all six panels. The colour coding in each 
panel represents a different galaxy property as follows: (a) Mass of the host halo (log scale), which is correlated with the environment 
density, (b) Morphology, quantified by the bulge-to-total mass ratio, (c) Mass-weighted stellar age. (d) Mass-weighted stellar metallicity. 
(e) Total stellar mass (log scale), (f) Star formation per year per unit stellar mass (log scale). Black corresponds to galaxies with no star 
formation at all. 
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Figure 11. Morphology variations in a model with no bulge cri- 
terion for cooling shutdown. The similarity to the bulge panel 
in Fig. 1101 indicates that the bulge criterion is of secondary im- 
portance to the correlations between colour, luminosity and mor- 
phology. 



Fig. El shows that this is indeed a transition from a central 
galaxy to a satellite. 

Fading into the red sequence need not be the end point. 
Red galaxies, in hot massive haloes, do not accrete gas, but 
they can still merge with other galaxies. The highest merger 
rates are for massive galaxies in the bright part of the red 
sequence, as higher mass implies stronger dynamical fric- 
tion. Thus the more massive a galaxy is when it leaves the 
blue sequence, the more it will continue growing by merging 
after it has joined the red sequence. Repeated mergers with 
other massive galaxies allow the formation of giant ellipticals 
with stellar masses of ~ lO 12 A/0, which would be unattain- 
able through simple accretion of gas in the presence of a 
limit on cooling. Most of this growth is due to dissipation- 
less merging within the red sequence. This interpretation 
of the role of merging and passive evolution is supported 
by the dominance of central galaxies in the bright part of 
the red sequence and of satellite galaxies in the faint side, 
separated roughly at L„, the characteristic luminosity of a 
galaxy at the top of the blue sequence (Fig. [3J . It is also 
supported by the variation of morphology with position in 
the colour-magnitude diagram (Fig. HUb t. Blue galaxies are 
mainly spirals, with bulge-to-disc ratio that increases with 
luminosity. The morphological composition along the red se- 
quence also varies from mostly spiral galaxies at M r ~ — 20 
to mostly elliptical galaxies at M r ~ —22.5, evidence for an 
increasing merger rate from faint to bright galaxies. 

The correlation of colour and magnitude with morphol- 
ogy is a generic outcome of the merger scenario and is not an 
artifact of the bulge criterion for shutdown. This is demon- 
strated in Fig. 1111 in which the bulge criterion was elimi- 
nated. The central galaxies of group and cluster size haloes, 
shown as red dots at the bright-end of the red sequence 
in Fig. [3] are those with the most intense merger history 
and therefore with the highest bulge-to-disc ratio. A quali- 
tatively similar correlation with morphology is valid already 
in the "standard" model without cooling shutdown. 

The presence of two populations within the red sequence 
manifests itself also in the environment panel of Fig. 1101 The 



galaxies that live in the most massive haloes (1O 13_14 M0) 
dominate both the bright end and the faint end of the red se- 
quence, but the middle range between M r = —21 and —22.5 
shows many middle-of-the-range haloes of 10 12_13 Mq and 
even smaller. This means that as the halo mass grows in 
time the galaxy can evolve according to one of the following 
two tracks: it may either keep on growing until it becomes a 
bright E/SO galaxy or fade to become an early- type satellite 
with a signifcant disc. This pre diction of our mode l is con- 
sistent with the SDSS results jBlanton et alJ [2006 ) , where 
the number density in the l/i -1 Mpc environment of galax- 
ies is measured to peak both at the bright end and at the 
faint end of the red sequence while the morphology shows a 
gradient along the red sequence. 

We learn from the star formation rate panel of Fig. 1101 
that model galaxies are red because they are not forming 
stars and thus lack a young stellar population. Indeed, the 
distinction between the red sequence and the blue sequence 
is predominantly due to age rather than metallicity, as can 
be deduced by comparing the gradients at a given luminosity 
in the age panel and the metallicity panel of Fig. 1101 How- 
ever, the colour gradients along the red sequence and along 
the blue sequence themselves are driven by the metallicity 
gradient more than by the age gradient. 



9 CONCLUSION 

We have addressed the origin of the robust division of the 
galaxy population into two major types: the blue sequence 
of gas-rich galaxies with a young stellar population and the 
red sequence of gas-poor galaxies dominated by old stars. 
While the blue sequence is truncated beyond a characteristic 
stellar mass of ~ 3 x 1O 1O M0, the red sequence extends to 
higher stellar masses. This bimodality is strongly correlated 
with the morphological type and with the galaxy density in 
the environment. Large galaxy surveys have revealed that 
the bimodality is reflected in almost every other property 
of galaxies and that the separation between the two types 
is fairly sharp. High-redshift surveys show that a fraction of 
the red sequence is in place already at z ~ 1 — 2, and that 
massive star-forming galaxies exist at z ~ 2 — 4. 

We find that models of galaxy formation can success- 
fully reproduce the bimodality once a complete shutdown 
of cooling and star formation is imposed in haloes above a 
critical mass of ~ 10 12 Mq starting at z 3. The time evo- 
lution of the bimodality is reproduced when allowing partial 
cooling and efficient star formation above the critical mass 
at higher redshifts. By incorporating these simple new fea- 
tures in hybrid semi-analytic/N-body simulations, we have 
achieved an unprecedented simultaneous match to the joint 
distributions of galaxy properties at low and high redshifts, 
practically confirming the qualitative predictions of Dekel & 
Birnboim (2006). 

The motivation for the sharp transition at a critical halo 
mass comes from theoretical analysis (BD03; lBinnevll2004l : 
DB06) and cosmological simulations (K05; DB06) showing 
a sharp transition from galaxy build-up by cold filamentary 
flows in small haloes to virial shock heating of the halo gas 
in massive haloes. The abrupt shutdown in the supply of 
cold gas for star formation, due to the abrupt appearance 
of stable virial shocks as opposed to the smooth variation 
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of cooling time with mass, is responsible for the very red 
colours of the red sequence, the colour gap between the red 
and the blue sequence, and the sharp truncation of the blue 
sequence above L*. 

An important input for the success of the model is that, 
once the gas is heated to the virial temperature, it never 
cools. In this model shock heating serves as the trigger for 
efficient AGN feedback, which maintains the gas hot. The 
long-term effect may be associated with a self-regulated ac- 
cretion cycle maintained by the coupling between the cen- 
tral black hole and the hot gas. The luminous quasar mode 
is associated with the rapid build-up of massive black holes 
probably due to the high gas inflow rate triggered by galaxy 
mergers. The resulting high output power may be important 
in termina ting star formation and ending the quasar phase 
itself (e.g- lDi Matteo et aljEoblh . but this phase is short- 
lived and cannot serve to maintain the gas hot for long peri- 
ods. The self-regulation is achieved when the accretion rate 
drops to an equilibrium value at which the energy ejected 
from the black hole and absorbed by the gas compensates 
for the radiative losses (Cattaneo 2006). The mechanical ef- 
ficiency of this mode is close to 100%, while in the quasar 
mode most of the output power goes into radiation. In the 
quiet mode, the dense clumps of cold gas are not swift in 
responding to the external injection of mechanical energy — 
they can be destroyed or blown away only when the AGN 
is sufficiently powerful as in a quasar. On the other hand, 
the hot gas, which feeds the black hole in a gradual fash- 
ion, is more dilute, and can therefore respond effectively to 
the injected energy. Therefore, the self-regulated mode of 
hot accretion may be responsible for a significant fraction 
of the mechanical energy affecting the gas, even if the hot 
gas contributes only a small fraction of the final black hole 
mass. 

In a parallel work, ICroton et al.l <|2006F) have imple- 
mented a related semi-analytic model, where they focus on 
the growth of black holes and AGNs as feedback sources. 
In particular, they assume a transition from rapid to slow 
cooling and associate the accretion of cold and hot gas with 
an 'optical' and a 'radio' mode of black-hole growth, respec- 
tively. The latter provides the relevant AGN feedback pro- 
cess and determines the gradual shutdown cooling in massive 
haloes. We focus instead on the thermal properties of the 
gas as the trigger for an abrupt shutdown above a threshold 
halo mass, which translates into the bimodality scale in the 
galaxy properties. The key to our successful modelling of the 
galaxy properties is this robust shutdown, regardless of the 
details of AGN feedback. This simple addition to the model 
results in an unprecedented detailed match to the data at 
low and high redshift. 

Our "new" picture introduces a strong correlation be- 
tween the spectral galaxy type and the host halo mass, while 
the morphological type is driven by the galaxy merger his- 
tory as in the "standard" model. Thus, while the new model 
gas fraction and colour of massive galaxies are different, the 
morphological distribution and the correlation between mor- 
phology and environment density are reproduced as in the 
standard model. The natural correlation between the mass of 
the host halo and the density of the enviro nment, according 
to the "halo model" o f galaxy distribution ( Yan et al]l2003l: 
lAbazaiian et al.ll2004l : iKravtsov et al 1 I2004F 7 is reflected in 
the correlation between spectral and morphological type. 



A few comments on the evolution in the colour- 
magnitude diagram, as summarised in Fig. ll2l ('see also Dekel 
& Birnboim 2006). Spiral galaxies grow discs by cold fila- 
mentary streams at the centres of haloes below the critical 
mass. As the accreted gas is converted into stars, the galax- 
ies become brighter and redder along the blue sequence. A 
galaxy leaves the blue sequence either when its halo grows 
above the critical mass or when it merges into such a halo 
and becomes a satellite galaxy. In both cases, the galaxy 
ceases to accrete gas from its hot environment and rapidly 
turns red and dead. Small satellite galaxies tend to fade 
as they redden and turn into the red discs observ ed at the 
faint end of the red sequence (IBlanton et al. 2006). Big cen- 
tral spirals whose haloes have grown above the critical mass 
continue their growth along the red sequence through dis- 
sipationless mergers and evolve into the giant ellipticals at 
the bright end of the red sequence. Gas-rich mergers may 
also exhaust the cold gas reservoir of the merging galaxies 
and transform two blue spirals into a red elliptical, but this 
is a less frequent path in the formation of massive early- 
type galaxies. This picture implies that the colour difference 
between the blue and red sequences is primarily due to stel- 
lar age, while the colour-magnitude gradient within the red 
sequence is largely a metallicity effect, since deeper poten- 
tial wells retain their metals more effectively. The colour- 
magnitude gradient along the blue sequence is a combination 
of the two effects. 

A detailed study of the origin of the red sequence via 
evolutionary tracks in the colour-magnitude diagram is de- 
scribed by Cattaneo, Dekel & Faber (in preparation). At 
high redshifts, the enhanced star-formation rate that we 
have incorporated in galaxies of all masses helps curing the 
difficulty of the standard scenario in forming enough stars 
in massive galaxies at z ~ 2 — 4, as indicated by the obser- 
vations of LBGs and SCUBA sources. This new ingredient 
is motivated by the build-up of galaxies via nearly super- 
sonic filamentary flows (BD03; K05; DB06). The collisions 
of these streams among themselves and with the central discs 
are likely to trigger efficient bursts of star formation, in anal- 
ogy to collisions of discs or cold clouds. However, the fate of 
cold streams in hot haloes and the resulting star formation 
are yet to be studied in detail. 

We note that our scenario helps explaining the "down- 
sizing" of the formation of massive galaxies, where more 
massive galaxies tend to form stars earlier and over shorter 
periods. The most massive galaxies, those above the shock- 
heating mass, have formed their stars at z > 2 — 3 by cold 
flows in hot media, and stopped forming stars at lower red- 
shifts. Galaxies below the critical mass continue to form 
stars at lower redshifts. 

While the simple prescriptions that we have used for 
star formation and its shutdown are clearly only crude ap- 
proximations, their success in simultaneously matching all 
the observed bimodality features at different redshifts indi- 
cates that they fairly represent the effects of the complex 
physical processes involved. This motivates a detailed study 
of the processes of cold-flow induced star formation and the 
subsequent shutdown above the critical mass at late times 
due to the coupling between virial shock heating and AGN 
feedback. 
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Figure 12. A schematic sketch of galaxy evolution. Galaxies grow 
along the blue sequence through cold filamentary flows until they 
stop accreting gas because their host halo has grown above the 
critical shock-heating mass. At that point they move into the red 
sequence by becoming red and fading in luminosity. The shutdown 
of cooling above the critical halo mass sets an upper limit to 
the luminosity of blue galaxies and explains the characteristic 
galaxy luminosity L*. Galaxies get to the E/SO bright end of 
the red sequence either by dissipationless ('dry') mergers along 
the red sequence, or by gas-rich ('wet') mergers from the top of 
the blue sequence. This simple cartoon ignores the complications 
associated with the temporary reddening by dust extinction in 
star-forming galaxies. 
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